ANNALES UNIVERSITATIS MARIAE CURIE-SKLODOWSKA 
LUBLIN-POLONIA 
VOL.XLIL1 SECTIO AAA 1987 


Zaktad Fizyki Technicznej 
Akademia Rolnicza w Lublinie 


Krzysztof MURAWSKI, Roman KOPER 


Model Equations for Waves Propagation in an Infinite Cylinder 
о with Thin Walls of Elastic Rings 
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1. INTRODUCTION 


Compliant tubes containing a fluid have been studied since the time of Thomas 
Young [1] in connection with modelling the propagation of the arterial pressure 
pulse. Moodie [2] et al. has discussed the tube equations based upon a thin- wal- 
led shell theory for thethered tubes. The fluid has been assumed inviscid and a 
one-dimensional. theory has been extracted by averaging quantities over the tube 
cross section. The other model has been employed to study the propagation of pres- 
sure and flow pulses along initially uniform tubes and their subsequent interaction 
with various junctions characteristic of the arterial system [3]. Wave propagation 
and shock formation in nonlinear elastic and viscoelastic fluid filled tubes for a 
Mooney-Rivilin material has been discussed [4]. Also radial motion of a non-linear 
viscoelastic tube has been studied [5j. Two-dimensional analysis was employed to 
study pulse propagation in thin-walled circularly cylindrical elastic tubes contai- 
ning an inviscid and incompressible liquid [6]. A viscoelastic shell theory model for 
transcient pressure perturbations in fluid filled tubes has been presented and tested 
against experiments involving water filled latex tube [7]. 

In 1980 Lamb [8] showed using Taniuti-Wei's [9] method that.in the absence 
of energy dissipation the fundamental set of equations describing liquid that is 
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confined within an infinitely long circular cylinder with thin walls of elastic rings 
leads to the Korteweg-de Vries equation. For this case the same equation has been 
derived via the Lagrangian method [10]. On the contrary when dissipation has 
been taken into consideration Burgers equation has been derived !11, 12]. 

In the present paper we combine the above mentioned equations for a weak 
dissipation and dispersion case obtaining the Korteweg-de Vries-Burgers equation. 
Also the Korteweg-de Vries equation is derived via the derivative expansion method 
[13] in section 5. For a strong dispersion case the nonlinear Schrodinger equation 
is obtained via the derivative expansion method in section 3 and via the reductive 
Taniuti- Wei's inethod in section 4. Last section is devoted to the short summary. 


2. FUNDAMENTAL SET OF EQUATIONS 


We consider now the one-dimensional irrotational fluid waves of characteristic 
amplitude l and characteristic length А in an infinitely long tube with thin walls of 
elastic rings and a diameter 2a to take into account dissipation of energy, nonline- 
arity and dispersion of medium on the assumption that | << 2a << A. Then the 
fundamental set of equations may be written as follows — , 

(i) the continuity equation 


(04), + (V A). = 0, (2.1) 
(ii) the Euler's equation 


УЕ + VV + ра = Vas (2.2) 


(ii) the Newton’s equation [9] 


E 2xa n( Eh — 2aq) 
—. = —- ео а 2.3 
"wt арт ^ Pmh” Pmh ( ) 
(iv) the state equation 
p= р(р) =4-р, (2.4) 


where the following notation are used: р — liquid density, А — area of the crossection, 
V - liquid velocity, м — viscosity coefficient, a — tube radius at the undisturbed 
uniform state, pm — density of the tube material, E - Young's modulus in the 
circumferential direction, p — liquid pressure, 4 — outside pressure, h — thickness of 
the wall. The subscripts z and t imply partial differentiation. 


3. THE NONLINEAR SCHRÓDINGER EQUATION. 
THE DERIVATIVE EXPANSION METHOD 


We derive now the nonlinear Schrodinger equation via the derivative expansion 
method [13]. For simplicity the effect of damping is neglected aud we assume that. 


i 
f 
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the fluid density is constant. In other context, Davey [14] has derived a nonlinear 


Schrodinger equation which is modified to allow dissipation. 
Introducing new coordinates as follows 


1 2a 2ap 
As 228 EL a V —\/—Vv, 
mur eas tak Eh 
E 2p 
t — ; | E. z= ——r, 3.1 
V рта? рт ah | ) 
dimensionless equations are obtained, 1.е., 
А, + (AV), 20, 1 (3.2) 
V, + УИ, + А, + Ав, =0. (3.3) 


We define two dimensionless small parameters, namely: 


e=, (3.4) 
l 
= 23 Я (3.5) 


Е and ó measure the weakness or dispersion and nonlinearity, respectively. The 
nonlinear Schrödinger equation is derived on the assumption that & = ô. 
In equations (3.2) - (3.3), we introduce the multiple spatial and temporal scales 


In ==", =", n*-0,1,2..., 


zo, to mti. (3.6) 
The derivative operators are considered to be of the form 
à д 8 20 
erat pee E rien соли 1а SiT 
Ot Ato tei Ae 8t4 чӣ s 
à à à ,8 E 
an ieee T mts. ee (3.76) 
929 Эдж ^ X Otoütiüzo 912942; / 
^ 93 93 93 93 
ва ee RR y Qc ЕР ыы SRI. 3.7 
Š (sias 900920 E Oa, 98924 )+ ort 


The dependent variables V, A are expanded around the undisturbed uniform state 
into the asymptotic series in terms of the same parameter 5 by writting 


А=1+6А; + 52А +:::, У = 6V, + 67Vo+---, (3.8) 
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Substitution of (3.7) and (3.8) into equations (3.2) and (3.3) yields a sequence of 
equations by equating the coefficients of like powers of Е: 


Alto ar Viz, =0 , (3-9) 
=: : 
Vie, + Ais, + Аз» =0, (3.10) 
Аз, + Art, + Voz, + Viz, + AiVizo + Vi Aia, 0, (3.11) 
ri 


Vot, + Vie, + УИ + А22, + Aiz, + Азаа, + 2А}. 0,20 + Ais, = 0 , (3.12) 


Aat, + А, + Аш, + Vaz, + Vos, + Viza + А (Иа + Viz,)+ 
+A2Viz, ap Vy (A2z, ap Aiz,) $ V2A1s, =0 , (3.13) 


Vo, + Aszo + A3t3zo + Vat, + Vit; + Vi(Vas, + Vizi) + V2Vizo + Аза, t 


А12, + 2A2totizo + Аза, + Atia, + 24160625 + 2Aitotiz: + Ад», = 0. (3.14) 
The linear equations (3.9) and (3.10) may be solved to get 


Ay = $(zi, М, .. Je? + cc. + ол (а, 45, +.) › (3.15) 
Us em peer + с.с. + Fini. sts), (3.16) 
9 = Ехо —wto, (3.17) 
k2 
. 28 
а = р, (3.18) 


where с.с. stands for the complex conjugate to the proceeding term and із introdu- 
ced for the reality of A, and Vi. $ is a complex function and o; and fj are real 
functions which describe interaction between slow mode and wave trains [15]. This 
problem, however, is not discussed in this paper. We assume that a; = £4 = 0. 

Substituting equations (3.15) and (3.16) into the second order equations (3.11) 
and (3,12) we obtain 


Aot, + Vos, + CT + Фа) e + с.с. + 2iwd2e7” + сс. = 0, (3.19) 


Vet, En Arzo a A222, sp (Z4. ar $s, 3E 20Еф,, = u^.) e'? nasus 


us 
M d + с.с. =0. (3.20) 
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This set of equations may be solved to get 


1 - ; , 
Аз = zG* t cc. + b(zi,..., t5... Je? + ec. Emus tn. ), (3.21) 


Vz = a a ам + с.с. + HET seesi...) te.c. + G(z1,...,¢1,-..) (3.22) 
where b, f ‚ 6, are arbitrary functions which satisfy the following relation 
i(wb — kf) = du + in (3.23) 
and $ satisfies the characteristic equation 


фе: + was, =0. (3.24) 


Finally, from the third order equations (3.13) and (3.14) the nonlinear Schrodinger 
equation may be obtained 


Р i 2u? (1 — k? 
2iw (1 + Е?) фе, + ze — k?) + К2]ф., – Or Pes 
LJ 
=k? фа + фа, + tw (1+ 22), + iw fe, + А, +b2,)+ (3.25) 
—iku?b,, = ш(ёш + 269)ф . 
If we introduce coordinate transformation defined by 
& = т — wets , Ty 72088) (3.26) 
and 
2 = 22—w&t2, т = ёз (3.27) 


from (3.24) we obtain that Kcu = 0 and all terms containing b and f are eliminated 


because of (3.23) and dispersion relation (3.18). Equation (3.25) may be written 
in the more convenient form 


ч 


и? (12k8 + 354 + 39k? 


Ф 7 ана SK? + 8) 


ee) =C, (8.28) 


mi 
2(1 + k?)? one: 
where C is an arbitrary constant to be determined from the boundary conditions. 

4. THE NONLINEAR SCHRODINGER EQUATION. 
TANIUTI-WEI'S METHOD j 


In the previous section, we have derived the nonlinear Schrodinger equation via the 
derivative expansion method. Here, turning our interest to the same case of strongly 
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dispersive waves, we develop similar analysis based on the reductive Taniuti-Wei’s 
method. For this aim we expand quantities A, V into the following series 


=1+ Le 5 Al? (e, refe) 


"m--—oo 

Pa + e^ NM vine, rje alkan wt (4.1) 

n=l m= — om 

where 

AQ) =v =0, mzzil, (4.2) 

Asn 
up = | ugh =u), з) 

vi 
ê= e(z — №), r=et. (4.4) 


The asterisk denotes the complex conjugate. Substitution of the expansion (4.1) 
into the fundamental set of equations (3.2), (3.3) yields the following equations for 
the n-th order terms 


[^] 

yi! = = AP | (4.5а) 

Е: 2 

k 

C Mu BE 
, 1+’ i 
-AAU) — iw A?) + V0 + ikv? =0, (4.6a) 
-AV(?) — ioV(?) + (1 — w? = 25) AU) + ik(1 — o?) AP?) = 0, (4.65) 
«2: 

| kV? — wA?) + w(Al)? = 0, (4.6c) 
- 209 VÍ? — 2k(4? — 1) A + = (ay =0, (4.64) 


AU! — ЛА — io Al?) + V (9 + ikv АКА" + (473) 
4.7a 


"io AUS AC) + if АЙ) + ky? AU) =0, 


(2 + awk) AU) — Av (2) — wl?) + io AC VO + уу ufu 


(1 — 02 — 29k) AP) + ik(1 — 2) APP (20A + KA?) AUS, = 0, 
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e: 
АФ! = аис и +3) ано ie. 
ee в 
From equations (4.6a) - (4.64), we obtain 
T: -Л А - so AP! tiv P — 0, (4.8) 
; A= We, (4.9) 
Al?) = E (al)? , (4.10) 


yp Saas А 2 (4.11) 


A compatibility condition for the components of ys and AÜ' is reduced to the 
nonlinear Schrödinger equation for the first order quantity 


fone! MAE a WU NIS gla 
"bo den capi р ee О 


2 
= spa (e AD? + 2kV(?)) ACD. | (4.12) 
Substituting (4.7c) and (4.7d) into (4.12), we obtain 


w(12k6 + 35k* + 39k? + 9) За T 
salt) 4 1 AQ 2 AQ AY = AY, 
4k4(k4 + 3k? + 3) Mi Aa 2(1+ k2)2 tee 
we 
where Ó is a new constant defined by 
GE Cla es + 2kC3). (4.14) 


2k? 
Cı and Cz may be calculated from the boundary conditions. 


5. THE KORTEWEG-DE VRIES EQUATION 
DERIVED VIA A DERIVATIVE EXPANSION METHOD 


Our main purpose now is to apply a derivative expansion method in order to 
obtain the Korteweg-de Vries equation which describes the propagation of small 
amplitude and long waves. We introduce the multiple spatial and temporal scales 


bs SCO; zu Е т. (5.1) 
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The dependent variables V, A are expanded around the undisturbed uniform 


state by writing ` 
oo oo 
Va S P WV. © 4А=1+ 5" An: . (5.2) 
n=1 n=l 
The Korteweg-de Vries equation is derived on the assumption that 
==”. к (5.3) 


Substituting (5.1) and (5.2) into the fundamental set of equations (3.2), (3.3), we 
get a sequence of equations by equating the coefficients of like powers of e. The 
first three sets of equations may be written in the followig form 


Ain T Viz. =0, (5.4a) 

=? ; 
Vi, T Ais; == 0 , (5.45) 
Ait, + Viza =), (5.5а) 

e»: = 

Vit, SP Aiz, =0, (5.55) 
A21 ap Ait; gis Van, F (41Vi)z, + Vis; =0 1 (5.6a) 

=“: 
Vas, + Иль, + ViViz, + Ада, + Ais, + АЙ =O. ‘ (5.65) 

From equations (5.4) and (5.5), we find 

У, = Vi(&1 = 21 — tı) = Ai(&), (5.7) 
Vi = Vi (f2 = 22 — t2) . (5.8) 


The fourth equations (5.6) lead to the following equation 


Vas 7 Va + Artsa, t (Aiii) + (И Ае Је + Vias (5.9) 
Vias F(Viin)a + 41660 + Ares 70 


The second order terms may be removed if we assume that V2 depends on z, and 
t; through £1. Thus we have 


2415 + 2Aiz, + 3A; Are, + Алё, ёё = 0. (5.10) 


Transforming to the coordinate system moving with the phase velocity equals to 1, 
1.е., 
Ёз = тз із, т=з, (5.11) 


we obtain from (5.10) the Korteweg-de Vries equation 


3 1 
Air + 541416 + 5416566 70- (5.12) 
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6. THE KORTEWEG-DE VRIES-BURGERS EQUATION 


Our aim here is to obtain an approximate single model equation which describes 
the behaviour of small amplitude and long waves on assumption 6 = =2. For this 
purpose we apply a nonlinear perturbation Taniuti-Wei’s method the idea of which 
is to introduce the following coordinate stretching 


€ = Ve (z — Vot), reet. (6.1) 


Since we consider weakly nonlinear waves, we expand dependent variables around 
the undisturbed uniform state as power series in terms of the same parameter є: 


9+ єр + e) pa +..., 
Ao +=А, +24: +., (6.2) 
eV; + eV +... " 


<ь 3 
нии 


where Ао = ха2. 
Here we make an assumption that viscosity coefficient is small and proportional 
to the parameter \/Е 
в = УП . z (6.3) 


For this case, thus we consider a balance between the nonlinearity, dispersion and 
dissipation of energy. 

Introducing (6.1) — (6.3) into the fundamental set of equations (2.1) — (2.4) and 
equating all the coefficients of the various powers of є, we obtain for the first power 
of e: 


1 1 
A, = nhhT гац = $, (6.4) 
S 
Vimus ME л 6.5 
° 4(45 + Ао51) e 
where 
E 2ra $1 (952 AP AoS1)Vo 
$,2 ——, S=-—, =>) ==. “ices 6.6 
> Ohm Dials he RS 2405; 3% 
Finally, from =2, we obtain the Korteweg-de Vries- Burgers equation for ф: 
br + 8% + Мфее + apee =0, (6.7) 


where the nonlinear £, the dissipative M and the dispersive а coefficients are defined 


as follows: . 

= (Аоф + 2Vo) (992 — Уя) (6.8) 

992 Ао + (9 + 91 4o)Vo 

PFA 
M= m o. —- 6.9 
4[(ч+ gı Ao)Vo + 992 Ao] i 
4 

az Что (6.10) 


51 [(9 + 91 Ао) У + 992 Ао] 


Ш 
[2] 
[s] 
[4] 
"Bl 
fe 
Ч 
[8] 
[9] 
[10] 
ш] 
[22] 


[13] 
[14] 
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7. SUMMARY 


Basing on the derivative expansion and the reductive Taniuti-Wei's methods, 
it was shown that the fundamental set of equations is reduced to the nonlinear 
Schrodinger equation or the Korteweg-de Vries and the Korteweg-de Vries - 
Burgers equations depending on if the system is strongly or weakly dispersive. The 
Korteweg-de Vries-Burgers equation is derived on the assumption that viscosity 
coefficient is small and proportional to the small parameter =. 

In a surface water waves context the nonlinear Schródinger equation has'been 
derived and studied in some details [16]. It is well known that this equation can be 
solved with the help of the inverse scattering transform [9], possess Painlevé test 
and Backlund transformation [17] to obtain an N-envelope soliton solution and 
periodic envelope one [18]. 

It is worth noticing that the two coefficients by the second and third terms in 
equations (3.28) and (4.13) are responsible for modulational stability of waves [16]. 


If their product is negative nonlinear wave solutions are modulationally stable. In 


our case it occurs for arbitrary value of k. So, waves in tubes are modulationally 
stable. 

The Korteweg-de Vries- Burgers equation and its solutions has been discussed 
in some detail, see e.g. [19] and [20]. Also the Korteweg-de Vries equation has been 
modified to include both the dissipative and dispersive effects of viscous boundary 
layers [21]. In the frame of the two-dimensional theory of long gravity waves a mo- 
del equation was derived which combines the Kadomtsev-Petviasvili with Burgers 
equation [22]. 

The Korteweg-de Vries equation (5.12) agrees in form with the same equation 
obtained via other methods [9, 10]. 
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STRESZCZENIE 
W pracy zastosowano teorie fal nieliniowych, opartą na metodzie redukcji Taniuti- Wei i me- 
todzie wielu skal dla wyprowadzenia nieliniowego równania Schrodingera, równania Kortewega- 
de Vriesa і równania Kortewega-de Vriesa-Burgersa. Równania opisują propagacje fal w płynie 
wypelniajacym cylinder o cienkich i elastycznych pierscieniach. 


РЕЗЮМЕ 


В работе применено теорию нелинейных волн, основанную на методах pe- 


дукции Танюти-Веи и методе многих параметров, для получения нелинейного 
уравнения Шрёдингера, уравнения Кортевега-Де Фриса и уравнения Кортевега- 
Де Фриса-Бургерса. Эти уравнения описывают распространение волн в жидкости, 
заподняющей цилиндр с тонкими и упругими кольцами. 


Złożone 26.X1.1986 


